Skip to content

Encoding of female mating dynamics by a hypothalamic line attractor

Liu, Nair, Coria et al. (2024)

DisciplineNatural science
MethodExperiment
Tagscalcium imagingdata · experimentfemale sexual behaviorhypothalamusline attractormethod · machine learningmethod · regressionmethod · simulationneural coding

Summary

Problem: 雌性动物的交配接受性(receptivity)通常被视为由脊柱前凸(lordosis)反射定义的二元行为,但实际上雌性接受性在交配过程中及不同生理状态下均表现出高度动态性。然而,与雄性相比,雌性对成功交配动态的重要贡献长期被低估和研究不足。此前研究虽已鉴定出控制雌性接受性的神经回路(如VMHvl中的Esr1+神经元),但大脑如何编码交配过程中雌性行为的动态变化仍基本未知。

Approach: 作者使用微型显微镜(miniscope)对自由行为的雌性小鼠VMHvl中Esr1+、Npy2r−(称为“α细胞”)的钙活动进行纵向成像,结合无监督动力学系统建模(rSLDS)、广义线性模型(GLM)、光遗传学扰动、行为学视频注释等多种手段,系统分析了交配过程中雌性行为的动态特征及其神经编码机制。

Finding: 研究发现VMHvl α细胞群体活动在交配过程中呈现一个近似的“线吸引子”(line attractor)动力学结构。该吸引子包含一个具有长时程常数(约110秒)的“整合维度”,其活动在交配过程中持续斜坡式上升,在射精前达到峰值。该维度的活动水平与雌性接受行为程度正相关。光遗传学扰动实验证实了吸引子的真实存在——群体活动被推离吸引子后会弛豫回到吸引子上。纵向成像显示,该吸引子动力学随动情周期出现和消失,并依赖激素状态。

Significance: 该研究首次在自由行为的哺乳动物中证明了下丘脑线吸引子编码一种持续的、渐增的雌性性唤起或驱力状态,为理解大脑如何编码动态内部状态提供了新的理论框架。同时,该研究还首次证明吸引子动力学可被激素状态在数天时间尺度上可逆调节,拓展了人们对神经吸引子可塑性的认识。

Theoretical Framework

N/A

Research Design

本研究为实验性研究,采用多种互补的实验设计:

  • 行为学分析:对雌性小鼠与雄性交配互动的视频进行人工标注,识别十种雌性运动行为,并分类为食欲性(appetitive)、接受性(accepting)和抵抗性(resisting)行为。
  • 钙成像实验:使用微型头戴式显微镜对自由行为的雌性小鼠VMHvl中Esr1+、Npy2r− α细胞进行纵向钙成像,记录交配过程中的神经活动。
  • 动力学建模:使用无监督循环切换线性动力学系统(rSLDS)对神经群体活动进行建模,识别低维动力学结构。
  • 光遗传学扰动:对VMHvl网络进行光遗传学抑制,同时进行钙成像,检验吸引子的“吸引力”性质。
  • 行为扰动:通过光遗传学激活雄性VMH中的Sf1+细胞诱导雄性防御状态,人为中断交配,检验雌性吸引子活动的持续性。

分析单位:单个神经元(单细胞水平)和神经元群体(群体水平)。 关键自变量:雄性交配行为(嗅闻、爬跨、插入)、雌性交配行为(接受、抵抗、食欲性行为)、动情周期阶段、激素状态。 关键因变量:VMHvl α细胞的钙活动(z-scored ΔF/F)、整合维度的活动水平、线吸引子结构的存在与否。 操作化方式:行为通过视频人工标注分类;神经活动通过miniscope钙成像记录;吸引子通过rSLDS模型的时间常数谱和线吸引子评分识别;接受性程度通过接受行为(脊柱前凸和扭动)的时间比例量化。

Data & Sample

  • 样本量:行为学分析n=74次交配试验,n=28只雌性小鼠;钙成像n=15只雌性小鼠(每只16–207个神经元,平均89±12个);光遗传学扰动n=3只小鼠;纵向成像n=4只小鼠。
  • 动物:实验室雌性小鼠(具体品系未在摘要中明确,但使用了Esr1flp、Npy2rcre等转基因品系)。
  • 地区:美国加州理工学院(California Institute of Technology)。
  • 数据收集方法:视频行为标注、微型显微镜钙成像、光遗传学结合成像。
  • 时间周期:交配互动观察5–10分钟;纵向成像覆盖动情周期多个阶段。
  • 数据来源:原始实验数据。

Analytical Strategy

  • 行为学分析:计算各行为类型的概率随时间变化、行为持续时间分布、行为转换概率等。
  • 单细胞分析:使用选择概率(choice probability)评估单细胞对特定行为的选择性;使用广义线性模型(GLM)预测神经活动,评估行为变量和细胞耦合对神经活动的解释力(交叉验证R²)。
  • 群体动力学分析:使用rSLDS(无监督循环切换线性动力学系统)拟合神经群体活动,计算各维度的时间常数,识别“整合维度”;计算“线吸引子评分”(两个最长时间常数的log2比值)。
  • 验证分析:使用监督式目标降维方法(targeted dimensionality reduction)和FORCE训练的循环神经网络独立验证吸引子结构。
  • 扰动分析:光遗传学抑制期间观察神经轨迹偏离吸引子及恢复情况;使用未扰动数据拟合的模型预测扰动后轨迹。
  • 行为解码:训练线性解码器区分交配回合与回合间期(IBI)的视频帧。
  • 相关性分析:整合维度活动与接受行为比例的相关性。
  • 对照条件:雄性被限制在带孔容器中(pencil cup)以排除持续雄性气味对持久活动的贡献。

Results & Findings

1. 雌性交配行为是动态且分级的,而非二元 雌性交配行为并非简单的二元接受/拒绝,而是呈现动态变化:接受行为概率逐渐增加,抵抗行为先增后减。雌性83.9%的交配行为时间用于响应雄性(male-responsive),仅16.1%为自发行为(self-initiated)。雄性自发交配行为时间是雌性的11.3倍,表明雄性主导交配互动。在交配回合间期(IBI),雌性仍持续表现出社交行为,提示存在一种持久的内部接受/参与状态。

2. 单细胞水平上,VMHvl α细胞对特定交配行为的选择性较弱 仅约2–13%的α细胞对特定交配行为表现出“调谐”(choice probability >0.7或<0.3),大多数细胞表现出“混合选择性”。GLM分析显示,仅用雌性行为只能解释约8%的神经活动方差,加入雄性行为后约15%。然而,加入细胞耦合项后,解释方差大幅提升至约46%,表明局部回路相互作用对神经活动方差的贡献大于行为本身。

3. 群体水平上存在一个线吸引子,包含一个长时程“整合维度” rSLDS分析发现一个时间常数显著长于其他维度(110.7±13.6秒 vs 24.5±5.1秒)的“整合维度”,形成线吸引子结构。该维度在交配过程中持续斜坡式上升,从雄性嗅闻、爬跨或插入开始时启动,跨越多个交配回合和IBI持续上升。约30%的单细胞表现出超过25秒的自相关半宽(ACHW),提示持久活动。整合维度的斜坡活动是群体属性而非单个神经元的简单叠加——不同细胞在不同时间达到峰值。

4. 线吸引子具有真实的“吸引力”性质 光遗传学抑制VMHvl网络期间,神经轨迹被推离吸引子(偏离流形),抑制结束后轨迹回到吸引子上最近的固定点。模型对未扰动数据的预测能准确预测扰动后的轨迹。在不同吸引子位置施加扰动均观察到恢复现象,证实了吸引子上不同固定点的吸引力。

5. 整合维度活动在交配中断期间保持稳定 在自然交配回合间期(IBI)中,整合维度活动保持在与交配回合期间相似的高水平,无法用该维度训练解码器区分交配回合与IBI。在人为诱导的长时间交配中断(1–5分钟)期间,整合维度活动仍保持升高数分钟,表明该状态的维持不依赖持续的雄性接触输入。

6. 整合维度的斜坡活动与接受性程度正相关,且不单纯反映机械感觉输入 整合维度活动与接受行为比例正相关。在雄性插入率高但雌性接受行为低的试验中,整合维度几乎没有斜坡活动,表明斜坡活动并非简单地累积来自雄性插入的机械感觉输入。

7. 吸引子动力学随动情周期和激素状态可逆出现和消失 纵向成像显示,线吸引子动力学在动情周期的接受期出现,在非接受期消失,且依赖激素状态。这证明吸引子可被激素在数天时间尺度上可逆调节。

8. 意外发现 与先前使用批量钙成像的研究(显示VMHvlCckar神经元活动从爬跨开始到射精单向下降)不同,本研究发现整合维度在交配过程中持续上升,提示群体水平的动力学可能揭示出单细胞或批量记录无法捕捉的编码机制。

Limitations

根据论文文本,作者承认或暗示的局限包括:

  • 行为学分析依赖人工视频标注,可能受主观判断影响。
  • 钙成像记录的是α细胞亚群(Esr1+、Npy2r−),不能完全代表VMHvl所有神经元的活动。
  • 光遗传学扰动实验样本量较小(n=3只小鼠)。
  • 线吸引子是“近似的”(approximate),并非数学上严格的线吸引子。
  • 研究主要在小鼠中进行,跨物种推广性有待验证。
  • 论文文本中未完整列出所有局限性(部分内容在正文截断处之后)。

Key Contributions

  • 首次在自由行为的哺乳动物中发现并验证了下丘脑VMHvl中的线吸引子动力学,证明其编码雌性交配过程中的持续、渐增的内部状态。
  • 提供了线吸引子在自由行为动物中的直接扰动证据(光遗传学抑制结合成像),此前此类扰动尚未在自由行为动物中实现。
  • 揭示了雌性交配接受性并非二元行为,而是由神经群体动力学编码的分级、动态状态。
  • 证明了神经吸引子可被激素状态在数天时间尺度上可逆调节,拓展了对吸引子可塑性的理解。
  • 整合了行为学、钙成像、动力学建模和光遗传学等多种方法,建立了从行为到神经动力学的完整因果链条。
  • 发现局部回路相互作用对VMHvl α细胞活动方差的贡献远大于行为变量本身,提示内部动力学在编码中的核心作用。

Key Claims

"Females exhibit complex, dynamic behaviours during mating with variable sexual receptivity depending on hormonal status. However, how their brains encode the dynamics of mating and receptivity remains largely unknown." (Abstract)

"These observations suggest that a hypothalamic line attractor encodes a persistent, escalating state of female sexual arousal or drive during mating. They also demonstrate that attractors can be reversibly modulated by hormonal status, on a timescale of days." (Abstract)

"The continuous 'ramping up' in the integration dimension observed over a long timescale during mating is unexpected; it contrasts with studies of bulk calcium activity in female VMHvlCckar neurons, which have shown a unidirectional decrease from the start of mounting until ejaculation." (Results, "Neural dynamics in female VMHvl" section)

"Definite evidence for the attractive nature of the fixed points discovered by the rSLDS requires performing neural perturbations orthogonal to the line attractor. Such perturbations for line attractors have yet to be performed for freely behaving animals." (Results, "Perturbations of the female line attractor" section)

"Together, these data indicated that VMHvl α-cell activity displays line-attractor dynamics during mating, and that although male contact is integrated along the line attractor, the stability of the system in the integration dimension does not require continuous male contact-dependent sensory input." (Results, "Perturbations of the female line attractor" section)


My Notes